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^ ■ Abstract 

o, 

CJ , We report the value of the dynamical critical exponent z for the six dimensional 
Ising spin glass, measured in three different ways: from the behavior of the energy 

Cn ■ and the susceptibility with the Monte Carlo time and by studying the overlap-overlap 

^—N , correlation function as a function of the space and time. All three results are in a 

cn [ very good agreement with the Mean Field prediction z = 4. Finally we have studied 

^^ ' numerically the remanent magnetization in 6 and 8 dimensions and we have compared 

f^ . it with the behavior observed in the SK model, that we have computed analytically. 
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1 Introduction 

Nowadays the only way to calculate analytically critical exponents in non exactly soluble 
models is the use of field theory techniques. One can choose to use the e-expansion [U| or 
work at fixed dimension P|. 

These approaches have had a great success in the study of the pure Ising model (i.e. 
the 0^ theory). 

In the spin glass case the starting theory is a 0^ theory in the limit of zero components 
of the fields (induced by the replica trick). Another theories that can be described by the 
0^ field theoretical description are the percolation, the Potts models and the Lee- Yang 
singularities. The upper critical dimension of this theory is 6. 

In order to check the validity of the predictions of the e-expansion (where the calcula- 
tions are easier than work at finite dimension) we need to have numerical results on the 
critical exponents in dimensions near the upper critical one [du = 6), i.e. to simulate in 5 
and 4 dimensions (where e = 1, 2). 

The previous discussion refers to the static of the system and can be extended in a easy 
way to the dynamics of the spin glasses. 

Recently two of the authors (G. Parisi and P. Ranieri)^] has been able to compute the 
one loop correction to the dynamical critical exponent z whose Mean Field (MF) value 
(the base of this e-expansion) is 4. They found 

.i^)=4(l-^), (1) 

where e = 6 — d. 

To test this analytical result is mandatory to do numerical simulations in five dimen- 
sions. 

This paper is the first step to the calculation of the dynamical critical exponent in five 
dimensions. We can found in the literature results in four [Q and three dimensions [^. 

In this paper we use different techniques in the six dimensional case (where we have as 
reference the Mean Field value) to extract the value of z. We plain in the future simulate 
in five dimensions using the methodology developed and tested in the present paper. 

In particular we obtain the z exponent using three different off-equilibrium ways: 
the decay of the energy, the growth of the non linear susceptibility and by means to the 
scaling of the overlap-overlap correlation function. These three methods provide us of three 
estimates of z or ratios of z with the static critical exponents in a very good agreement 
with the Mean Field predictions (z/mf = 1/2, rj^p = and Zmf = 4). 

Whereas the z calculations have been done at the critical temperature of the system 
we have performed extra numerical simulation inside the cold phase to monitorize the 
"expected" dependence on the temperature of the z exponent (like has been obtained in 3 
and 4 dimensions 0, ^) and to check the predictions of De Dominicis et al. (monitorizing 
the growth of the susceptibility) that predict that the propagator restricted to the g = 
ergodic component goes like p~^, where p is the momenta of the propagator. Obviously at 
the critical point we expect the usual dependence on the momenta, i.e. p~^. 

We have also studied an interesting observable, the remanent magnetization. This 
magnetization is defined as follows: we put the system under a large magnetic field, turn 



off it and follow the decay of the magnetization of the system. This observable has a great 
importance because its measurement is accessible experimentally 0]. This decay is like 

M(t) ~ r\ (2) 

which define the A exponent 0. Following the results of D. Fisher and H. Sompolinsky |jlO[ 



we expect that A will be equal to the MF value only for d > 8, i.e. for exponents which 
are reminiscent of the "magnetic field" the upper critical dimension is eight, not six. 

We check this fact numerically and analytically in the present paper obtaining the value 
of the A exponent in six, eight and infinite dimensions (the SK model). 

2 Numerical simulation and observables 

We have simulated the ±1 six dimensional Ising spin glass, whose Hamiltonian defined in 
a hypercubic of volume L^ is 

where < i,j > denotes sum to nearest neighbor pairs, Jij = ±1 (with the same probability) 
are quenched variables and Si = ±1 are spin variables. 

We have simulated L = 8 and L = 10 sizes at the critical point with a number of 
samples respectively of 103 and 13, and also 106 samples of an asymmetric system of size 
12 X 8^ 

The static of this model was studied by Wang and Young [|ll[|. Simulating lattice sizes 
up to L = 8 they found that the static critical exponents were compatibles with the Mean 
Field predictions (z/mf = 1/2, ''7mf = 0) and that there were logarithmic corrections to the 
Mean Field exponents, effect linked to the upper critical dimension. Their estimate for the 
critical temperature was Tc = 3.035 ± 0.01. 

The main aim of this paper is to measure the dynamical critical exponent z in order 
to compare it with the Mean Field results (^mf = 4). To do this we have measured the 
behavior of the energy and susceptibility as a function of the Monte Carlo time 

E{t) - Eoo oc t-^ (4) 

X{t) oc t" (5) 

and the q — q correlation function. 

The first observable that we will examine is the dependence of the energy with the 
Monte Carlo time. We assume that at the critical point (and only at the critical point) 
it is possible to connect the approach to equilibrium of the energy and of the equal time 
correlation functions to the equilibrium static and dynamical exponents. For example in 
the case of the energy we find: 

E(t) - Eoo oc r^^-(^)/^ T = r„ (6) 

where z is the dynamical critical exponent, dim(ii^) = d — \jv is the dimension of the 
energy operator and d is the dimensions of the space. Assuming d = 6 and v = z/^f = 1/2 
we have that the exponent of the energy decay at T = T^ is 6 = A/z. 



Analogously we can write for the non linear susceptibility 

x{t) = L'WW), Q{t) = Y^T.^^itMt), (7) 

i 

where a and r are two real replicas with the same quenched disorder, the following depen- 
dence on the Monte Carlo time 

x(t) « t(2-'')/^ r = T„ (8) 

as t ^ Tcq{L), where Teq(L) is the equilibration time, which should diverge as Teq(L) oc L^. 
Where we have used that dim(x) = 2 — 7]. Assuming the Mean Field value for rj we have 
that the exponent of the susceptibility is h = 2/z at T = Tc. 

From these formulae it is possible to calculate the dynamical exponent via two different 
observables. In the six dimensional case, ii z = Zmf = 4 we should expect a behavior like 
t~^ for the energy and like t^'^ for the non linear susceptibility. 

A third way to calculate the dynamical exponent is using the overlap-overlap correlation 
length at distance x and time t defined by 



G{x, t) = —J2 Wi+xn+x(^iri) , (9) 

i 

where again a and r are two real replicas with the same disorder. In the simulations we 
start from a random configuration (T = oo) and suddenly we quench the system to Tc 
or below. Then the system begins to correlate itself and we can define a time dependent 
off-equilibrium correlation length, ^(T, t), as the typical distance over which the system 
have already developed correlations different from zero, i.e. G{x,t) ~ for x > ^(T, t). 
The growth of this correlation length with the Monte Carlo time defines the dynamical 
exponent z trough 

^{T,t)oct'/<^\ (10) 

We have seen that in three and four dimensions 0, Q the data fit very well the following 
functional form 

G,M)^^exp{-(^)]. (11) 

Thereby, this will be the third way to obtain the dynamical critical exponent. This third 
estimate of z is independent of the values of the static critical exponents. 

The equilibrium overlap-overlap correlation function constraint to g = was obtained 
by De Dominicis et al. |]12[, which in 6 dimensions reads 




CsRSB(a;)|g=o ~ S ^-2 \iT <Tc ^^'^^ 



The fact that the equilibrium correlation function C{x) has a power law decay also at 
T < Tc, implies that spin glasses are always critical in the glassy phase (i.e. below Tc) 
and therefore we can relate the off-equilibrium behavior of the correlation function to the 



equilibrium critical exponents. Using that the susceptibility is the integral of the correlation 
function 

X = Jd'xCix) (13) 

and that in the region where the susceptibility grows, following a power law of the time, 
the overlap is very small, we can use equations (|12D in equation ([13D , obtaining 

r 4I/2 if T = T 

If we take the limit T ^ Tc in the above equation we obtain that h{T), the exponent of 
the growth of the susceptibility, must be discontinuous at the critical point (i.e. h{T~) = 1 
while h{T+) = 1/2). 

Moreover, if we assume that l/z{T) is proportional to the temperature (this happens 
in 3 and 4 dimensions [0, 3) we must obtain a linear dependence on the temperature for 
h{T) in the low temperature phase. 

Finally we will study the decay of the remanent magnetization defined as 

1 ^' 
M(t,t^) = — ^(Ti(t)ai(t^), t>t^. (15) 

^ i=i 

This observable decays like 

M{t,t^) (xt-^. (16) 

We find that the Mean Field prediction for this exponent is A = 5/4. 

3 Results in d = 6 

The great part of the simulation work have been done at the critical temperature, chosen as 
the weighted mean between the one found by Wang and Young [|ll| (Tc = 3.035 ±0.01) and 
the one calculated by the series expansion [|13] (T^ = 3.027 ± 0.005): f3c = 0.3302 ± 0.0005. 
In particular we have tested that the exponents we measure do not vary in a significant 
way if the temperature is changed by an amount of same order than the uncertainty on 
Tc. To this purpose we have simulated the same system (of volume 8^) at the inverse 
temperatures Pi = 0.330 and P2 = 0.331, checking that the dynamics were compatible. 

Verified that, for the range of time and sizes we have used, the exponents we are 
interested in do not depend on the precise Tc choice, we have decided to run all the 
subsequent simulations at /3c = 0.330. At the critical point we have simulated more than 
one hundred samples of size 8^ and 13 samples of size 10^. The last number of samples 
may appear too small to average over the disorder; in fact we have used mainly the data 
from the 8^ systems to calculate the moments of the distribution of the overlaps. The data 
from the bigger systems have been used to measure almost-self-averaging quantities like 
energy whose fluctuations are very small considering that we are working with system with 
a million of spins. 

The results obtained are shown in figure |1| for the energy decay and in figure |^ for the 
non linear susceptibility growth. 
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Figure 1: Energy decay in a 10^ system at T = T^. 

We have tried to fit the energy decay both with a power law {E{t) = E^o + At~^) and a 
logarithmic law {E(t) = Eoo + A[\n(t/r)]^^). We are interested in the asymptotic behavior 
of the decay; then we fit the data in the range t G [tmin, oo) for various choices of tmin 
and we expect that the parameters of the fit converge rapidly when we increase tmin- The 
impossibility of fitting all the data with a single law (for t < 6) is due to the existence 
of an initial short time regime of a few steps during which the dynamic does not follow 
yet the asymptotic behavior. We find that the logarithmic law do not describe well the 
data because, even if it has more adjustable parameters, the best values of the parameters 
depend strongly on tmin, they are very correlated and they tend to unphysical values. On 
the other hand we find that fitting with the power law the values of the parameters Eoo, 
A and 6 converge to a stable value, with tmin of order of few MCS. In figure |l| we plot the 
power fit obtained with tmin = 6 (this is the lowest value for which the fit holds the x^ test); 
the best parameters are: E^ = -1.8880 ± 0.0001, A = 0.37 ± 0.01 and S = -0.98 ± 0.01. 
We note that the decay exponent is compatible with the mean field value {6mf = — !)• 

The line plotted in figure ^ is the best power fit to the susceptibility data. We have to 
take care when we try to interpolate in such a way because we know that the susceptibility 
growth follows a power law only in a limited time window; in fact at the beginning of the 
simulation the dynamics needs some short time to reach the asymptotic regime f\ and then 
because of the finite size of the system they have to converge to some finite value, i.e. the 
data of figure |^ have to converge to a plateau. These two effects may induce systematic 
deviations in the estimation of the h exponent: x(^) oc t^. In our case the first transient 
is almost absent thanks to the sufficiently high temperature, but problems may arise from 

^This initial time increases when the temperature is lowered. 
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Figure 2: Susceptibility growth in in systems of size 8^ at T 
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the existence of a plateau. In fact we find that fitting all the data points with the power 
law the x^ value is not acceptable while, just discarding the last few points, the fit works 
very well. The line reported in figure ^ is the power law obtained fitting the data in the 
range t G [6, 2000] 0; we note that the last three points are below the best fit, which is 
probably due to the effect just described. Our estimation of the dynamical exponent h is 
h = 0.49 ± 0.02 which is compatible with the mean field value (/imf = 1/2). 

We also show in figure |^ the results (all obtained on a 8^ x 12 lattice and with 200 
samples) for the h exponent in the low temperature phase {T/Tc = 0.5, 0.625, 0.75, 0.875). 
We also plot in this figure the value that we have obtain at T^ {h = 0.49). It is clear 
that h{T) is a discontinuous function at the critical point and that the limit from below, 
assuming a hnear behaviour, {h{T~) ~ 0.9) is almost twice the value of h(T^) = 0.5, in 
very good agreement with the correlation functions found by De Dominicis et al. [|12|. The 
quite small discrepancy can be due to the crossover between the two regimes in a finite 
lattice or to logarithmic corrections. 

Moreover we can see that the dependence of this exponent is well described by a linear 
law of the temperature, according to a z{T) inversely proportional to T. 

Measured the equal time q — q correlation functions defined in eq.(|^), we have studied 
the spatial correlations by the following technique (already used with success for the data 
analysis in 3 and 4 dimensions 0, §]). We expect a functional formula for this correlation 
function of the form 



x^ 



XiT,t)^ 



^Thc fits done in different time windows give similar results. 
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Figure 3: /i(T) exponent vs. T in the region T < T^. 

where, as usual, ^(T, t) oc t^/^C^) is the dynamical correlation length. For each value of the 
distance we fit the data of the q — q correlation function by the formula 



(j(a;,t) = 6*00(3^) exp y4(a;)t 



I--B 



V fixed X, 



(18) 



verifying that the value of the B parameter is almost independent from x and then fixing it 
during the following study. The dynamical exponent can be expressed as the ratio z = '~f/B, 
where 7 may be estimated by the power law fit A{x) oc x'^. By this way also our third 
estimation z = 4.2 ± 0.2 is compatible with the Mean Field value. In figure ^ we have 
plotted lnG(a; = 2,t) versus t, together with the best fit (eq. (|T8|) ). 

In the infinite time limit the function G{x,t) converges to Goq{x) which must give 
information on the q — q correlation function at zero overlap, calculated for the SK model 
by De Dominicis et al. |]T2|. They found that Gq=Q{x) oc x'^ at the critical temperature. 

By our simulation we find the GA^) plotted in figure |^ together with the best power 
fit in the range x G [2, 5]. We note that the point in a; = 1 if far away from the fit because 
for ln(a;) -^ — cxo the data must converge to 1 (by definition G{x = 0) = 1) while the power 
fit diverges. Anyway we are interested in the asymptotic behavior which seems to be well 
described by a power law, Goo{x) oc x~" with an exponent a = 4.2 ± 0.1, in agreement 
with the Mean Field result previously cited. 

4 Remanent magnetization 

The last part of this study has been dedicated to the decay of the remanent magnetization. 
We prepare the system with all the spins up (M(t = 0) = 1) and then we let it evolve 
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Figure 4: lnG'(x = 2, t) vs. t for the asymmetric system at T = Tc 
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Figure 5: Goo{x) for the partially asymmetric system at T = Tc 



toward the equilibrium where, in absence of any external field, no magnetization have to 
remain. In the cold phase we expect the decay of the remanent magnetization to be of 
algebraic type (see equation (|D); in particular we are interested in the exponent of such a 
decay at the critical temperature (which hereafter will be called simply A), to compare it 
with the same exponent of the SK model. 

We report in the following subsections our calculation of the A exponent in the MF 



approximation (see reference |T^ for another calculation of the A exponent in the spherical 
spin-glass model), together with the numerical verification and the estimations of such an 
exponent in the finite dimensional cases {d = 6 and d = 8) and also for the SK model. 

4.1 Analytical and numerical results in the SK model 

To study analytically the dynamical properties of the SK model we define the following 
soft-spin Hamiltonian: 

m=-PY: J^,S^S, + 1 ^ ,2 ^ 1^ ^ 4, (19) 

ij i ■ i 

where Sj {i = 1 , . . . , A^) , are one-dimensional real variables and Jij is a symmetric matrix 
with independent elements following a Gaussian distribution with zero mean and variance 
proportional to 1/A^. From the random matrix theory |T^ we know that, in the thermo- 



dynamic limit (i.e. A^ goes to infinity), the probability distribution for the eigenvalues of 
Jij is given by the semi-circle law: 

a(/i) = ^(4-^2)i/2 1^1 < 2. (20) 

A relaxation dynamics is introduced by the Langevin equation: 

where ^j(t) is a Gaussian noise with zero mean and variance {^i(t)^j(t')) = 26ij6(t — t'). To 
study the dynamical evolution of this model, we diagonalize the Jij matrix and we consider 
the dynamics of the projections s"'(t) of the spin Sj(t) on the eigenvector directions ip^ (with 
eigenvalues //„), such as Sj(t) = J2n^"(t)xl'^, where n = l,--,N is the eigenvector index. 



The properties of independence and orthonormality of the eigenvectors, |jT5| let us define 
the following Langevin equation for the component s": 

dts-it) = (/?/.„ - i)s"(t) - 1 E s°(t)/(t)s^(t) Ecv'f^z • c + at). (22) 

'^- Q/37 i 



In the thermodynamic limit, as shown in IM, the Hartree-Fock approximation is exact 
and we can linearize the Langevin equation: 

dts^it) = (/?/x„ - l)s-{t) - |c(t, t)s"{t) + Cit) (23) 
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where C{t,t) is the dynamical auto-correlation, 



C{t,t') = {s,{t)s,{t')), (24) 



evaluated for t = t'. We represent the average over the thermal noise as ((■ ■ ■)), while (■ ■ ■) 
indicates the average over the disorder as usual. 

With respect to the eigenvalues the auto-correlation function can be defined as follows 



C{t,t')=J df,a{fi){s^{t)s^{t')). (25) 

The self-consistent solution of the eq. ( PBD is 

gn^^^-^ ^ _gn(^Q^|g(/3M„-l)tg-(;/2/Jdt'C(t^t') _^ j* ^f ^{lii^u-l)(t-t") ^-9/2 f^,, dt' C{t\t') ^n^^lf^ ^ (^26) 



ft 
/O 

where t = is the initial time. 



We want to analyze the evolution of the system at T = Tc from a uniform initial 
condition: s"(0) = 1 , Vn. From ( PU] ) we obtain the following self-consistent equation for 
C(t,t): 



C(t, t) = H\t) J d^ a(/z)e^(^-^)* + 2 /^ dt" (^-^ j / d^ a(/z)e^(^--^)(-*") (27) 

where 

i7(t) = e-^/2/o^*'^(*'*'). (28) 

Let us suppose for H(t), at T = T^ a time-dependent asymptotic behavior (t -^ cxd) like: 

H{t) ~ tPe-^\ (29) 

where p can be determined auto-consistently from eq . (|27|) . This hypothesis implies for 
C{t,t), at large t, the following behavior 

C(t,t)~(2 + AT,)-^. (30) 

where a is an opportune constant. The critical temperature of the Hamiltonian (|T9|) is 
given by Tc = T° -|- AT^ = 2 + AT^ (T° = 2 is the critical temperature of (19) when g = 0) 



and that a perturbative calculation gives ATc = —2g and so 2/3c = 1 + g [|I^ . 

We recall that for T r^ T^. the bigger contribution to the dynamical relaxation of the 
spins comes from the region of the maximum eigenvalue, n = 2. 

At the critical point, the first term on the right side of the eq. (pTf ) scales with a power 
law t^^/'^+'^p. To be consistent with the hypothesis (|29| ) and (30) it should be p < 1/4. 



Now, we have to estimate the second term, which is proportional to: 

2j['dt"/d/.a(/.)^e-^/^^(^-'^)(-*") (31) 

If we consider t' = t — t" and we define the exponential as integral in the complex plane, 

^-2/3.(2-,). =Y^t^ [2/3,(2 - /.)tf = /^ 1^ n-s) [2/3.(2 - /.)tT , (32) 

11 
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Figure 6: Integration path. 



where C is the path shown in figure |^ (i.e. s = sq + ir, where r E TZ, and sq whatever real 
number in (—1,0)) and r(s) is the Euler Gamma function, we left with: 



2/^ fdt'Ti-s)- 
Jc 2m JO 



t^P 



-t" 



{t - t'fp 

After the integration over t' and /i we obtain: 
8 r ds 



d/ia(/i)(2/3J^(2-/i)^ 



(33) 



-v/tt Jc 2Tii 



T -s)t 



i+,r(s + i)r(i-2p)^,^2;5j^r(s 



3/2) 



r(s-2p + 2) 



r(3 + s) 



(34) 



To evaluate the integral ( ^4]) we analytically continue the function on the left of the path 
C, i.e. in the region where Re s < sq. Thereby, we have to consider the residues of the 
poles in this region. The residue of the pole in s = — 1 gives the constant contribution 
to the auto-correlation function C{t,t) while the time-dependent behavior, for large t, is 
determined by the value of p. In fact, for p = 1/4, r(s + 3/2) is simplified by r(s — 2p + 2), 
and, for large t, the leading behavior of C{t, t) comes from the residue of the pole at s = —2: 



C{t,t) ~ const +t 



-1 



(35) 
-3/2 and we 



consistently with the hypothesis ( p9D and (|30|). 

For p 7^ 1/4, on the contrary, we do not manage to cancel the pole at s 
should obtain 

C(t,t) -const + t"^/^ (36) 

in contrast with the previous hypothesis. Thus the only consistent solution for H(t) ~ f'e"* 
is p = 1/4. 
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An other way we can obtain the value of the p exponent is by solving the eq. (27) in 
the Laplace transform. In terms of the function g{t) 



git) 



-2gt 



n(t) ■ 



(37) 



the eq. (^) at T = T^ becomes: 



',.„ s(t") 



«W^/d,.(„e-»<-...2^d<"^/ 



d/icr(/i)e" 



-2/3,(2-/i)(t-t") 



(3^ 



By using the asymptotic form (^) of C(t,t), we obtain the following asymptotic equation 
for the Laplace transform of g{t), that we will denote g{s): 



~9{s) 



+ 



d/icr(yu) 
dficr{fi) 



1 




s + 2/?,(2 - 
1 


-/^). 


s + 2/?,(2 - 


-/i)J 





(39) 


ATA ~, . CL f°° , ~, n1 
2 j^(^) + 2A ^^^(^)J 


(40) 



By averaging over the eigenvalue distribution we obtain: 



9[s) 



^2s\ I 

+ 



2(3, A0V V2/3c 4/3, 



3/2 



AT, 



A 7 \ n r^ 



X] 



(41) 

By remembering that (3, = 1/(2 + AT^), we can developed the previous equation in 
ATc. We will also assume that in the limit s ^ 0; J^ dxg{x) is negligible with respect to 
g{s), and we finally obtain 

~g{s)r^^ + 0(1). (42) 

thus, for t -^ oo, g(t) ~ l/t^'"^ and from (|37D we obtain p = 1/4. This solution implies 
that J^dxg{x) < g{s). 

At this point, we can determine the decay rate of the correlation between the system 
and the initial conditions, i.e. the scaling decay of the remanent magnetization: 



M{t) = C{t,0) = {s,{t))= di2a{p){s^{t)) 



d/xa(/i)/J(t)e(^^-i)* 



t3/2 



^1/4 



-5/4 



(43) 
(44) 



Then, the analytical prediction for the exponent A, such as M(t) ~ t^^, is, in mean field 
limit, A = 5/4. 

The result is quite different from the one obtained in the spherical model at T = Tc 
where A = 3/4 0. 



For a numerical confirmation of this result we have simulated three SK models at 
criticality (T^ = 1) of sizes A^ = 480, 992, 2016 with number of samples of 10000, 5000 and 
1000 respectively, obtaining three estimations of the Amf exponent all compatible with the 
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Figure 7: Remanent magnetization in the SK model at T = T^ 



theoretical prediction. Since the data for the remanent magnetization have non evident 
finite-size effects, we have plotted if figure |^ the data averaged over all the simulated 
samples. The observable we have measured is M(t, t^ = 3) defined in eg. (PTSD , which follows 
the same decay of M{t) but has some advantages as will explain below. The line in figure |^ 
is the best power fit which gives Amf = 1.22 ± 0.02. 

We have repeated this numerical calculation for a temperature below the critical one 
(T = 0.8 Tc) and sizes N = 480,992,2016,4064. In this case we expect that the remanent 
magnetization tends to a non- zero asymptotic value due to the finite size of the system. 
So we have fitted the data of M{t, t^ = 3) with the following formula 



m(t, A^) = moo (N) + A r^^^^ 



(45) 



where we let A to depend on the temperature. Via a preliminar three parameter fit we 
have estimated A = 0.785(10) and with no systematic dependence on the lattice size. Then 
fixed the value of A to the just found value, we extrapolated the value of moo{N) by using 
a more simple linear fit, like the one plotted in figure ^ 

Using the values of moo{N) found by the previous analysis we were able to fit them 
by a power law of the system size: moo{N) oc N~'^''^^^^' . The data with the best fit are 
reported in a double- log scale in figure |^ (see ref. |T^ for a detailed study). 

Assuming a linear dependence of the exponent A(T) with the temperature, which has 
been observed in |]1^ [1^ for the SK model and in |T^ for a spin glass system on quenched 
0^ graphs (which should behaves like a mean-field SK model), and from the fact that 
A(T = 0.8) ~ 0.8, we obtain that the A critical exponent as a function of the temperature 
is discontinuous at the critical point (i.e. A(T~) ~ 1 while A(T+) = 5/4). 
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Figure 8: Remanent magnetization in the SK model at T = 0.8 Tc- 
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Figure 9: moo(iV) vs. N in the SK model at T = 0.8 T^,, 
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Figure 10: Remanent magnetization in ci = 6 at T = T^. 

4.2 Numerical estimation in d = 6 and d = 8 

The measurement of the decay rate of the remanent magnetization can be more difficult 
than one can imagine, because of the following effect: we know that M(t = 0) = 1, but we 
try to fit the M{t) data with a power law which diverges at t = 0. This effect is evident in 
a log-log scale where a power fit behaves like a straight line, while the M{t) data tends to 
the value 1 when Int — ;> — oo; in such a situation we have to discard the first data points 
to be sure we are measuring the asymptotic behavior. Unfortunately the useless data are 
just the ones with smaller relative error, while the data fitted are affected by a greater 
statistical indetermination which makes worst the estimation of A. 

One possible way to overcome this source of error is to measure some other observable 
that has the same behavior of M(t), but with a signal stronger and less fluctuating. Starting 
with all the spins up, the magnetization at time t is nothing more than the overlap between 
the conflguration at time t and the initial one. If we measure the overlap between the 
configuration at time t and the one at a fixed small time (t^ = 3 in our case), we expect 
that M{t,tu]) behaves like M{t), with similar statistical fiuctuations, but with a signal ten 
times greater. 

The results of the simulations of the EA model in 6 dimensions can be found in figure p!0| 
where we have plotted M{t, t^ = 3) versus the simulation time t; the line is the best power 
fit which gives an exponent A = 0.995 ± 0.005. This value is compatible with 1, but not 
with the Mean Field value Amf = 5/4. 

As explained in the introduction we believe that the upper critical dimension becomes 
du = 8 for such quantities linked, in some sense, to a magnetic field. In this case we start 
the simulation with the system totally magnetized, like if it was feeling an infinite strength 
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magnetic field, so the remanent magnetization may be one of sucli quantities. 

Tfien we fiave done some simulations of the ±1 Ising Spin Glass model in rf = 8 to see 
if we recover the Mean Field behavior of the remanent magnetization. 

The critical temperature in d = 8 has been extrapolated from the critical temperatures 



in d=3 [^1,4 [^ and 6 |[TTI|. In the limit oi d —^ oo the critical temperature diverges like 



Tc{d) ~ \'2d. Moreover, in the Bethe-Peierls approximation, there is an exact formula for 



the critical temperature |22, 23 



{2d-l)tanh\pc) = l (46) 

We have used this formula, which turn to be valid in the d -^ oo limit, as a starting point 
adding to it some terms which may mime the finite dimensions corrections. In particular 
we have substitute the r.h.s. of eq.(^) with a fourth order polynomial in 1/d (the term 
of zeroth order being always 1) and we have tried to fit the known critical temperatures 
fixing two term of the polynomial to zero and letting free the coefficients of the other two 
terms. Among all the 6 possible choices we have selected the one with smaller value of x^? 
which reads 

(2d - 1) tanh2(/5,) = 1 + ^ + ^ (47) 

with B = 0.95±0.14 and D = 117±4. Through this fit we estimate the critical temperature 
for the eight-dimensional EA model as f5c{d = 8) = 0.270 ± 0.001. The error reported is 
an underestimation of the real one because there are systematic deviations due to the 
arbitrarily choice of the fitting function. 

We have also repeated the analysis looking at the quantity T^ /{2d — 1), which takes 
values in the range [0, 1]. Knowing that in the d ^ oo limit holds that 

rp 2 

1 , (48) 



2d-l 

we have tried to fit the critical temperatures adding to the r.h.s. of eq.(^8D a polynomium 
in 1/d, obtaining as the best resulting interpolation the one reported in figure |TT], which 
reads 

^ = l-5,16(26)-L-5(l)J3. (49) 

From the plotted fit we obtain an estimation of (3c{d = 8) = 0.270 ± 0.001. Also this 
error should be increased because systematic deviations due to the arbitrarily choice of the 
fitting function are not taken into account. 

Anyway the estimations of (3c{d = 8) obtained with different interpolations are all in 
the range [0.260, 0.270] and we think that the true critical temperature is with very good 



probability in this range; in fact, as you can see in figure |TI|, the value of (3c{d = 8) is 
strongly dependent on the value of l3c{d = 6), which is known with high accuracy, and on 
the way l3c{d) tends to zero as the dimensionality is increased. 

For example if we assume that the successive improvements of the Bethe-Peierls approx- 
imation tend to increase the value of l3c{d) for each d, then the formula (|46|) will give a lower 



bound for the inverse critical temperature; in d = 8 this lower bound is Pdd = 8) > 0.264. 
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Figure 11: Polynomial fit to the critical temperatures. 
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Figure 12: Remanent magnetization in d = 8 at T ~ Tg. 
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From figure |TT] we can get also another significant information: the point where the 
fitting function crosses the x axis may give us an estimation of the lower critical dimension, 
which turn to be di ~ 2.65. 



In figure |T2| we have plotted the data, with the best power fits, of the remanent mag- 
netization in c? = 8 at temperature f3 = 0.260 and (3 = 0.270. From the fits we get the 
exponents A(/? = 0.260) = 1.256 ± 0.08 and X{/3 = 0.270) = 1.235 ± 0.013, which are both 
compatible with the MF result (Amf = 5/4). 

5 Conclusions 

We have calculated numerically for the first time the dynamical critical exponents in six 
dimensions in three different ways, all compatibles within the statistical error. 

Thanks to the previous results we can check also the static critical exponents (for 
instance getting the z value obtained from the scaling of the q — q correlation function), 
and we obtain values that agree very well with the static critical exponents and the critical 
temperature found in the literature |]ll] . 



Finally we have calculated analytically in the Mean Field approximation the exponent 
of the remanent magnetization and we obtain numerical results that confirm that for this 
observable the upper critical dimension is eight and not six. 

We plain in the future extend this work on the five dimensional Ising spin glass. 
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